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Employing daily wind data from the ECMWF, we perform passive particle ad- 
vection to estimate the Lagrangian velocity correlation functions (LVCF) associated 
with the midlatitude tropospheric flow. In particular we decompose the velocity field 
into time mean and transient (or eddy) components to better understand the nature 
of the LVCF's. A closely related quantity, the absolute dispersion (AD) is also ex- 
amined. 

Given the anisotropy of the flow, meridional and zonal characteristics are con- 
sidered separately. The zonal LVCF is seen to be non-exponential. In fact, for 
intermediate timescales it can either be interpreted as a power law of the form t~ q 
with < a < 1 or as the sum of exponentials with differing timescales - both in- 
terpretations being equivalent. More importantly the long time correlations in the 
zonal flow result in a superdiffusive zonal AD regime. On the other hand, the merid- 
ional LVCF decays rapidly to zero. Before approaching zero the meridional LVCF 
shows a region of negative correlation - a consequence of the presence of planetary 
scale Rossby waves. As a result the meridional AD, apart from showing the classical 
asymptotic ballistic and diffusive regimes, displays transient subdiffusive behaviour. 



I. INTRODUCTION 

Being a quantity of fundamental interest in the statistical characterization of a turbulent flow, 
the Lagrangian velocity correlation function (LVCF) has been the subject of detailed study [1], [2]. 
Indeed, if we view the midlatitude free troposphere to be in a non- ideal turbulent state [3], [4], it 
is of natural interest inquire into the nature its LVCF. Similarly, apart from being a measure of 
transport and its close relation to the LVCF [5], [2], the absolute dispersion (AD) is known to yield 
information about the structure of the underlying flow [6] , [7] . 

Ofcourse there exist other measures, such as the relative dispersion (or Lyapunov exponents for 
smooth flows) and finite scale analyses (or finite size Lyapunov exponents), which are a part of 
the complete statistical characterization of a flow. In particular, finite scale statistics are argued 
to be useful when the available range of scales is small and there is the possibility of crossover ef- 
fects from differing regimes [8] . Our approach is slightly different from earlier modelling studies 
of large scale atmospheric dispersion which employ a Langevin equation, equivalently assuming an 
exponential LVCF, to parameterize the dispersion process [9], [10] (a discussion of such an approach 
can be found in [11]). In the present note we will only consider the LVCF and AD, focussing their 
characterization and complementary nature. Also, we will see the accordance of these statistics with 
simpler physically relevant systems. 

For computational purposes, daily pressure level global T63 (192 x 96, 17 levels) resolution wind 
data from the ECMWF is used as the 3D advecting velocity field. We computed the trajectories 



X I thank Referee A for the discussion of finite scale statistics. 
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(and record the velocity along these trajectories) of an ensemble of passive particles. The trajec- 
tories were computed in latitude, longitude and pressure coordinates by a standard fourth order 
Runge-Kutta scheme and the velocity fields were interpolated in a linear fashion. Averaging is done 
with respect to an ensemble (denoted by S) of passive particles which remain in the free troposphere, 
i.e. particles which escape into the boundary layer or into the stratosphere are excluded from the 
statistics 2 . 



II. LAGRANGIAN VELOCITY CORRELATION FUNCTIONS 

A. Zonal LVCF 

Denoting the zonal velocity by u(A, <f>,p, t), the zonal LVCF (R u (t)) is defined as 



_ < u(x(t + r)) u(x(t)) > s 
Hu{ ) < u(x(t)) 2 >S 



Here x(t) represents the trajectory of an individual member of S. To account for the inhomogene- 
ity of the flow and to focus our attention the midlatitudes, S is chosen to comprise of N mem- 
bers located randomly such that they initially satisfy 0° < A(S) < 360°, 35° < <p(S) < 55° and 
400mb < p(S) < 700mb, where \,<j),p represent the longitude, latitude and pressure coordinates. 
Note that no restriction is placed on the trajectories, i.e. the statistics presented are not conditional 
on the particles remaining in the midlatitudes for all times. 

As can be seen in the upper panel of Fig. 1, which shows log[i? u (r)] for the winter and summer 
seasons, the zonal LVCF is clearly non-exponential. This is in agreement with studies on 2D [12] 
and geostrophic [13] turbulence, but stands in contrast to investigations of 3D turbulence [14], [15]. 
Further, for intermediate values of r (i.e. 2 < r < 10 days), it appears that (lower panel of Fig. 
1) R u {t) ~ T~ a with < a < 1 (a = 0.32 and 0.45 for DJF and JJA respectively). The issue of 
whether the correlation function is indeed a power law is known to be a delicate matter. Indeed any 
monotonically continuous function, such as the aforementioned power law, can be expressed (over 
a given range of scales) as the discrete sum of exponentials - see Berglund [16] for definitions and 
details - a fact which has been elegantly utilized in the analysis of queues in communication networks 
[17], [18]. Whether the zonal LVCF is the sum of two (or more) Ornstein-Uhlenbeck processes with 
differing timescales - as is argued to be an effective parameterization of the LVCF in 2D turbulence 
[12] - or, if it indeed shows power-law scaling is practically undecidable especially if the behaviour 
is seen over a small range of scales 3 . On the other hand, what is certain is the enhancement of 
the Lagrangian correlation time [6] (T u — J Q °° R u {t)dt). Indeed, as will be observed, the persistent 
correlations in the zonal velocity field result in a superdiffusive AD regime. 



Experiments were performed with an initial set of 20,000 and 30,000 particles, no discernable difference 
in the statistics was noticed in these two situations. 

3 Comments of Referee B are acknowledged as they provoked a more careful look at this issue. 
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B. Meridional LVCF 



Defined in a manner similar to Eq. (1), but using the meridional component of the velocity field 
(i>(A, <j>,p, t)) we compute the meridional LVCF (R v (t)). As is seen in Fig. 2, R v (t) decays to zero 
in one week. Interestingly, we notice a pronounced anticorrelation, i.e. R v (t) < before the final 
Rv(t) — > behaviour. On comparison with 2D [12] and geostrophic turbulence [13] we find this 
feature to be unique to the present situation. This is easily explained when one takes into account 
the latitudinal restriction imposed by the rotation of the planet [19]. Specifically, in the absence of 
strong diabatic or frictional effects, the conservation of potential vorticity gives rise to large scale 
Rossby waves [20] . We argue that the meridional oscillatory motion implied by the Rossby waves is 
responsible for the aforementioned anticorrelation. 



C. Eddy and Time Mean LVCF's 

To gain some insight into the connection between the nature of the LVCF's and the structure of 
the tropospheric flow, we partition the daily data into time mean and transient components. Specif- 
ically, the time mean is u(x) = (1/T) J a T u(x, t)dt (where T is the duration of the entire season) and 
the transient or eddy component is defined as u'(x, t) = u(x, t) — u(x). Even though the tropospheric 
flow does not posses a clear spectral gap, i.e. there is a near continuum of active (temporal and 
spatial) scales, we expect the above decomposition to separate processes which vary on scales that 
are farthest apart [21]. 

Fig. 3 and Fig. 4 show the zonal and meridional LVCF's computed by cxculsively utilizing the 
Eulerian time mean (upper panels) and Eulcrian eddy fields (lower panels). In spite of the crudeness 
of our partition, the similarity between the eddy LVCF's in both cases is evident. Indeed, apart 
from the slight anticorrelation retained in the meridional eddy LCVF, both R u >(t) and R v >(t) are 
rapidly decaying functions with a timescale of the order of a couple of days. 

On the other hand the zonal and meridional time mean LVCF's are strikingly different. Where 
Ru(t) is strongly correlated on long timescales - something we would expect from a slowly varying 
zonal jet - Rv(t) represents oscillatory motion as induced by a large scale wave. Moreover, com- 
paring the behaviour of R v < (r) and Ry (r) leads us primarily attribute the anticorrelation observed 
in R v (t) to the time mean component of the flow. Indeed, behaviour consistent with these results 
has been observed in studies of balloon trajectories in the Southern Hemisphere [22] (see especially 
their Figs. 9, 10 and the discussion regarding the timescales involved in the definition of stationary 
and transient components of the flow). 



III. ABSOLUTE DISPERSION 



The absolute dispersion (AD) is defined as, 

AD(i) =< (x(t) - x(0)) 2 >s (2) 
For ideal (i.e. isotropic, homogenous, stationary and zero time mean) flows, we have [5], [2], 

AD(t) = 4E [ R(t) {t-r) dr ; E = kinetic energy (3) 
Jo 

The short time limit (i.e. r — > 0) of the above yields ballistic motion whereas when T u is fi- 
nite, the long time limit (i.e. r >> T u ) yields diffusive behaviour [1]. The presence of other 
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exponents, i.e. AD(t) ~ V ; 7 7^ 1, is referred to as anomalous diffusion [6]. Before dis- 
playing the results, let us fix some notation. We denote the total AD by A(t), i.e. A(t) =< 
(x(t) — x(0)) 2 + (y(t) — y(0)) 2 + (z{t) — z(0)) 2 >$, where x,y,z are cartesian coordinates. Further 
the individual components of the AD are denoted by Ai{t) where i represents a coordinate, for eg. 
A x (t) =< (x(t)-x(0)) 2 > s . 

From Eq. (3) we have R(t) ~ t~" => AD(t) ~ t 2 ~ a . Even though the non-ideal nature of 
the present flow, especially its rich time mean structure [3], [4], is likely to invalidate the direct 
applicability of Eq. (3) - nonetheless, from the form of R u (t) (whether one takes it to be a power 
law or a sum of exponentials, both being equivalent from the discussion in the previous section) it 
is reasonable to expect the zonal AD to exhibit anomalous behaviour at intermediate timescales. 
Indeed as can be seen in Fig. 5 (which shows A x (t) as computed using the daily wind data) 4 , 

A x (t) ~ t 2 < t < 2 [Ballistic] (4) 
A x (t) ~ f ; 7 = 1.45 2 < t < 8 [Anomalous : Superdiffusive] (5) 

A x (t) ~ R\ ; R e = Earth Radius t>8 [Saturation] (6) 

In order to avoid the effective boundedness of the domain we unwrap the longitude and present 
A\(t) in Fig. 6. Now the anomalous regime (t s ; S = 1.6) lasts from Ti < t < T 2 days (Ti ~ 2 
and T 2 ~ 25), after which we see the beginning of an asymptotic diffusive regime (lower panel of 
Fig. 4) 5 . It is worth mentioning that the anomalous behaviour of A\(t) is in close agreement 
with laboratory experiments on quasi-geostrophic flows [23]. Further, zonal supcrdiffusion has been 
identified in studies involving large amplitude Rossby waves [24] and in more general PV conserving 
flows (where ofcourse, the superdiffusion is along PV contours) [25] 6 . 

Regarding the meridional AD, apart from the initial ballistic behaviour we expect to see normal 
diffusion at large t, as R v — > quite rapidly. Once again, we use Eq. (3) to get a feel for the 
meridional AD at intermediate timescales. Qualitatively approximating R v (t) ~ oT T l Cl cos(u;r) 
(Fig. 2), numerical integration of Eq. (3) yields the AD shown in Fig. 7. Apart from the two 
asymptotic regimes we notice transient subdiffusive scaling. This is in accord with results utilizing 
random shear flows, where anticorrelation in the LCVF was associated with subdiffusion and even 
complete trapping in extreme cases [26]. Indeed the actual meridional AD, A z (t) shown in Fig. 8, 



behaves in precisely the same manner, 

A z {t)~t 2 0<t<2 [Ballistic] (7) 

A z (t) - t 7 2 < t < 7 [Anomalous : Subdiffusive] (8) 

A z (t)~t t>7 [Diffusive] (9) 



4 A y (t) is virtually identical to A x (t). Also, A z (t) (shown later ) << A x (t) hence A(t) also behaves in the 
same fashion as A x (t). 

5 Note that we should expect S / 7 as A\(t) only involves changes in A whereas A x (t) is sensetive to both 
A and 4>. 

6 I thank Referee C for providing references on particle motion in PV conserving systems. 
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Comparing this behaviour with the meridional AD in strictly PV conserving flows [24], [25] we 
see that in those cases the meridional AD is bounded whereas the violation of PV conservation at 
large times (i.e. greater than a week) leads to unbounded normal diffusion in the present situation. 
Such superdiffusive zonal and subdiffusive meridional behaviour has also been recently observed in 
a study of anisotropic drift-wave turbulence [27] - particularly remarkable is the similarity of the 
anomalous exponents in the two situations. 



IV. SUMMARY 

Employing daily wind data from the ECMWF, we have estimated the zonal and meridional LVCF's 
of the midlatitude tropospheric flow. The zonal LVCF is seen to be non-exponential in character. 
Physically, given that the midlatitude tropospheric flow has a rich time mean structure along with 
an energetic eddy field [3], [4] - this observation is not entirely unexpected. Moreover, from this 
perspective our examination of Ru(t) and R u '{t) serves to clarify the roles of the time mean and 
eddy fields respectively. Specifically, the eddy field by itself generates an almost exponential rapidly 
decaying LCVF whereas the time mean component - roughly a slowly varying unidirectional jet flow 
[21] - is seen to be strongly correlated. 

Apart from decaying to zero in a relatively short time (w 1 week), the meridional LVCF exhibits 
an anticorrelation - i.e. R v (t) < before R v (t) — > 0. We attribute this anticorrelation to the pres- 
ence of large scale planetary waves - a basic consequence of PV conservation on a rotating planet. 
Examining R v >(t) we see that the meridional eddy LCVF is very similar to its zonal counterpart. 
Whereas Ro(t) - a manifestation of the large scale stationary waves - has an oscillatory character 
and indicates the time mean component to be primarily responsible for the above mentioned anti- 
correlation in R v (t). 

As regards the AD the point that stands out is the simultaneous existence of superdiffusive and sub- 
diffusive anomalous scaling in the zonal and meridional directions respectively. It must be stressed 
that we lack a quantitative relationship between the LCVF and the AD in this non-ideal situation. 
Nonetheless, a certain qualitative basis is provided by super- and sub-diffusive behaviour in ideal 
turbulent fields with enhanced (power laws or sums of exponentials depending on ones interpre- 
tation) and anticorrelated LVCF's respectively 7 . Finally, given that similar behaviour has been 
observed in drift wave turbulence [27], we are led to speculate on the possible universality of this 
phenomenon in fields where (slow) jets and waves co-exist with (fast) eddies. 
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7 Note that this anomalous behaviour is transient i.e., it is flanked on either side by an asymptotic regime. 
Even though in the present situation this behaviour is supported by the nature of the LVCF, it is worth 
keeping in mind that crossover effects (especially in the superdiffusive case) could play a role in determining 
the quantitative nature of the anomalous exponents [8] . 
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V. FIGURE CAPTIONS 



• Figure 1 : Upper Panel : Zonal LCVF. Lower Panel : Possible power law behaviour at inter- 
mediate timescales. Though a sum of exponential processes would result in similar behaviour 
at intermediate scales. 

• Figure 2 : Meridional LCVF : Note R v (t) < before R v (t) -> 0. 

• Figure 3 : Upper Panel : Time Mean Zonal LCVF (DJF data). Lower Panel : Eddy Zonal 
LCVF. Note the different timescales in the two panels. Also, by about one week the eddy 
correlations have almost completely died out wheras the mean flow is still strongly correlated. 
JJA data (not shown) behaves in a qualitatively similar manner. 

• Figure 4 : Upper Panel : Time Mean Meridional LCVF (DJF data). Lower Panel : Eddy 
Meridional LCVF. Once again, note the different timescales in the two panels. JJA data (not 
shown) behaves in a qualitatively similar manner. 

• Figure 5 : Zonal AD. Ballistic — > Superdiffusivc — > Saturation. 

• Figure 6 : Upper Panel : DJF Longitudinal AD (Ballistic — > Superdiffusivc — > Diffusive). 
Lower Panel : dA\(t)/dt Vs. t. Note dA\(t)/dt — > Const. => normal diffusion. 

• Figure 7 : Upper Panel : Synthetic LVCF R(t) ~ c- T / Cl cos(lot) (C\ =7,uj = 0.3). Lower 
Panel : Induced AD. 

• Figure 8 : Upper Panel : Meridional AD (JJA and DJF curves have been shifted for clarity). 
Ballistic -> Subdiffusive Diffusive. Lower Panel : A z (t)/t ~ t -> A z (t) / 't ~ t fi - 1 < (3 < 
-» A z (t)/t - Const. 
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FIG. 1. Upper Panel : Zonal LCVF. Lower Panel : Possible power law behaviour at intermediate 
timescales. Though a sum of exponential processes would result in similar behaviour at intermediate scales. 
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FIG. 2. Meridional LCVF : Note R v (t) < before R v (t) -> 0. 




FIG. 3. Upper Panel : Time Mean Zonal LCVF (DJF data). Lower Panel : Eddy Zonal LCVF. Note the 
different timescales in the two panels. Also, by about one week the eddy correlations have almost completely 
died out wheras the mean flow is still strongly correlated. JJA data (not shown) behaves in a qualitatively 
similar manner. 
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FIG. 4. Upper Panel : Time Mean Meridional LCVF (DJF data). Lower Panel : Eddy Meridional LCVF. 
Once again, note the different timescales in the two panels. JJA data (not shown) behaves in a qualitatively 
similar manner. 



11 




12 





13 




14 




15 



